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1.0 INTRODUCTION 


This report describes results from the Nozzle Thermal Analysis (NTA) contract 
NAS8-3961 1 which has been performed for NASA/Marshall Space Right Center in 
support of the Advanced Solid Rocket Motor (ASRM) development. The primary 
emphasis of this program has been to investigate regions of the ASRM nozzle which have 
potential thermal or surface erosion problems. This study has been performed with state 
of the art analysis tools developed by Aerotherm Corporation with specific application to 
solid rocket motors. Some of the tools are computer codes which are used throughout the 
industry and have become standards. Other of the tools have recently been developed for 
NASA under the Solid Propulsion Integrity Program (SPIP). Several areas of the ASRM 
nozzle have received special attention during this study. These areas include: 

1. The submerged region of the nozzle which contains the flex seal. 

2. The nozzle entrance region which is subjected to a severe particle impact environment. 

3. The material interface region in the nozzle exit cone. This is the region where the exit 
cone material transits from standard density to low density carbon phenolic. 

4. The aft region of the nozzle exit cone which may be subjected to particle impacts. 

Study of these nozzle regions has been performed for both the subscale MNASA 
test motor and the full scale ASRM motor. The approach which has been taken in this 
program is to perform calculations for the MNASA motor and compare results to 
measured data. These comparisons have been made to calibrate the calculations and 
understand the phenomena occurring in this motor. With models developed and 
calibrated, calculations were then performed to evaluate the ASRM nozzle, especially the 
critical regions mentioned above. This report contains results from these calculations and 
describes the assumptions and analysis procedures which have been used. 

The report is organized into six sections; section 2 contains a description of the 
MNASA and ASRM motors, their propellant and chamber operating conditions, Section 3 
describes the flowfield calculations which have been performed for both of the nozzles. 
Section 4 describes the thermal and erosion calculations which have been performed. 
Section 5 presents conclusions, and Section 6 presents recommendations for aspects of 
the ASRM nozzle which need additional attention. 
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2.0 MOTOR DESCRIPTIONS 


Figure 1 presents the motor and nozzle configuration used for the MNASA 
ASRM 48-5 test which is described in reference 1. A more detailed schematic of the 
nozzle is presented in Figure 2. This test motor has a short aft dome section which 
separates the motor chamber from the nozzle. The nozzle is a subscale test item with 
many features designed to simulate the full scale ASRM nozzle. The throat diameter of 
this nozzle is approximately 10 inches. 

The ASRM 48-5 test was designed to evaluate a large variety of carbon phenolic 
materials. The materials of greatest interest for the present study are MX4996 standard 
density carbon phenolic and FM5939 low density carbon phenolic. This nozzle also 
contained a flexseal insulation specimen which consisted of polyisoprene shimmed with 
FM5055 carbon phenolic. This specimen was nonfunctional but was placed in the 
submerged region of the nozzle to simulate the expected environment for the full scale 
ASRM. 

Figures 3 and 4 are similar presentations for the ASRM motor and nozzle. The 
major differences in the two configurations (i.e., MNASA vs. ASRM) are: 

1. the position of the propellant grain relative to the nozzle 

2. the location of the flex seal, and 

3. the size of the two nozzles, the ASRM has a 54.488 inch throat diameter. 

There are also many other less significant differences. The ASRM utilizes only two 
carbon phenolic materials, the standard density MX4996 and the low density FM5939. 

The propellant formulation presented in Table 1 is used by both of these motors. 
This is an 88/19 (%solids/%aluminum) HTPB propellant. Figures 5 and 6 present the 
pressure traces for each of these motors. Figure 5 is the measured pressure for the ASRM 
48-5 firing and Figure 6 is the estimated chamber pressure for the ASRM motor. 


FORWARD DOME IGNITER CASE SPOOL ALTERNATE AFT DOME 



PROPELLANT CARTRIDGE #1 PROPELLANT CARTRIDGE #2 

Figure 1. MNASA Motor Configuration 
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Figure 2. MNASA Nozzle Configuration 
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Figure 3. ASRM Motor Configuration 
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Figure 4. ASRM Nozzle Configuration 




Table 1. Propellant Formulation 
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Figure 5. Chamber Pressure for MNASA 48-5 Motor Firing 
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Figure 6. Estimated Chamber Pressure for ASRM 
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3.0 NOZZLE FLOWFIELD ANALYSIS 


This section describes the flowfield analyses which have been conducted for both 
of the motors. This includes nozzle heating rate and particle flow analysis. This section 
will present the results from these calculations as well as describe the procedures which 
have been used. Calculations for the MNASA motor will be presented first followed by 
those for ASRM. Unless otherwise stated the analysis procedures for the two are 
identical. 

3.1 MNASA Nozzle Flowfield and Heating Rate Calculations 

Nozzle flowfield, heating rate, and surface thermochemistry calculations have been 
performed with the SPIP developed Nozzle AeroThermochemistry (NAT, Reference 2) 
computer code. This code is a coupled version of several other computer codes. The 
ACE module of the code is used to perform one dimensional isentropic flow solutions, 
surface thermochemistry solutions and provide inputs to the Navier Stokes solutions. The 
Navier Stokes solutions were performed with the KIVA module. These calculations all 
used the "Mollier" option in the code. With this option the code treats the propellant gas 
as a single species which has real gas properties which are determined by the ACE 
module. These properties are those for the gas in chemical equilibrium. This option was 
used because it greatly increases the efficiency of the KIVA calculations. This chemical 
equilibrium assumption has been found to yield acceptable heating rate results for solid 
rocket motors, i.e., chemical kinetic effects are insignificant. Table 2 presents the 
elemental composition used for the propellant gases. Some of the trace elements, e.g., 
iron have been ignored. This table also presents the average chamber pressure for the 
MNASA motor and the corresponding chamber temperature. 


Table 2. Propellant Gas Composition 


Element ' ' 

Hi Fraction (%) 

Hydrogen 

41.58 

Carbon 

9.55 

Nitrogen 

6.82 

Oxygen 

27.27 

Aluminum 

8.06 

Chlorine 

6.71 

P c = 45.4 atm 

T n = 6418° R 


Figure 7 presents the analysis grid which was used to perform the flowfield 
calculation for the MNASA nozzle. This flow solution includes the aft dome as well as 
the nozzle submerged region. The analysis grid has been refined in regions where largest 
flow gradients are expected to occur. Results from the one dimensional isentropic flow 
solutions (ACE) were used as initial conditions and inflow conditions for the KIVA Navier 
Stokes calculation. The inflow conditions, temperature, density and velocity were applied 
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to the left hand side of the analysis grid. An ambient pressure of 1 atmosphere was 
applied as the outflow condition at the right hand side of the grid. 



Figure 7. MNASA Nozzle Flow Analysis Grid 

Results from the flow calculations are shown in Figures 8 through 11. Figure 8 
presents flow velocities calculated for the nozzle. Figures 9 through 1 1 present pressure, 
temperature and Mach contours, respectively. Each figure identifies the values of the 
maximum and minimum contours as well as the increment. A, between each contour. The 
two dimensional nature of the flow is greatest just down stream of the nozzle throat. 

These flow calculations for the MNASA nozzle required 2100 cycles. Each cycle 
is a small incremental step in time. These steps in time are necessary because the KIVA 
module is a time accurate Navier Stokes solver. A calculation will require a greater 
number of cycles to reach steady state depending on how far the initial conditions are from 
the steady state solution For this computation, the flow in the submerged region did not 
completely reach steady state conditions. This was due in part to the various convergence 
criteria used in the calculation (The flow in the submerged region is very low velocity and 
tends to oscillate due to small variations in pressure. The pressure convergence criteria 
used for these calculations was 1.0E-04. A smaller criteria is probably needed to totally 
eliminate the observed flow oscillations. However, the magnitude of these oscillations is 
small enough that the calculated heating conditions should bound the actual heating 
environments.) 
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= 0.142 


tours 


The results from the Navier Stokes solution were used as boundary layer edge 
conditions and the convective heating rates were then determined with the MEIT module 
in the code. This module uses law of the wall assumptions and a momentum energy 
integral technique (MEIT) to determine convective transfer rates. The boundary layer 
edge conditions are taken as the flow conditions calculated for the cell adjacent to the 
wall. The boundary layer starting point is specified by the NAT user and for these 
calculations was taken as the point deepest into the submerged region. The edge 
conditions calculated with the KIVA module are shown in Figures 12 through 14. These 
figures present the edge velocity, pressure and temperature as a function of streamlength 
measured from the boundary layer starting point. A summary of these results is also 
presented in Table 3. This table presents edge conditions as a function of station, axial 
position and streamlength. The nozzle stations will be identified later in Figure 20. 

A second analysis grid was developed to more accurately calculate the flowfield in 
the submerged region of the nozzle as well as the throat entrance regions. This grid, 
which is shown in Figure 15, is referred to as the "cutoff' grid since it is terminated 
immediately aft of the nozzle throat. This grid is significantly refined compared to the grid 
presented in Figure 7. Calculations were performed with this grid for 3000 cycles. Even 
still, the fluctuating flows in the submerged region were not completely eliminated. 
Nevertheless, the flow did "steady out" significantly so that the flow conditions in this 
region could at least be bounded. Figure 16 presents velocities calculated in the 
submerged region of the nozzle for two points in time. The velocities in the submerged 
portion of the nozzle are on the order of 50 to 100 ft/sec. 
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Table 3. MNASA Boundary Layer Edge Conditions 


Station 

Axial 

Position (in) 

Streamlength 

(in) 

Ve 

(ft/sec) 

Te 

(deg-R) 

Pe 

(atm) 

4 

3.550 

1.81 

31.52 

6401.3 

44.90 

5 

2.477 

2.91 

65.99 

6404.2 

45.00 

6 

1.393 

4.03 

72.81 

6398.2 

44.70 

7 

0.365 

5.09 

89.03 

6394.2 

44.60 

8 

-0.701 

6.19 

107.2 

6390.3 

44.40 

9 

-1.923 

7.44 

182.0 

6402.6 

45.00 

10 

-2.863 

8.42 

272.5 

6409.4 

45.30 

11 

-3.815 

9.48 

459.4 

6415.8 

45.70 

12 

-4.232 

10.50 

1400.0 

6373.9 

45.70 

13 

-3.974 

11.10 

1305.0 

6339.2 

40.90 

14 

-3.520 

11.80 

2069.0 

6271.0 

37.50 

15 

-2.361 

13.30 

3112.6 

6140.1 

36.60 

16 

-1.681 

14.00 

3981.9 

5990.2 

29.70 

17 

-1.171 

14.55 

4369.0 

5902.4 

26.28 

18 

0.0 

15.65 

5313.0 

5689.8 

19.80 

19 

0.981 

16.70 

6381.2 

5143.9 

12.06 

20 

2.064 

18.15 

6714.9 

4857.4 

7.281 

21 

2.972 

19.17 

7098.1 

4721.9 

6.208 

22 

3.882 

20.17 

7389.8 

4592.1 

5.507 

23 

4.802 

21.18 

7582.6 

4538.9 

4.923 

24 

5.504 

21.96 

7706.7 

4488.4 

4.458 

25 

6.172 

22.66 

7810.8 

4444.2 

4.191 

26 

6.888 

23.43 

7917.4 

4398.2 

3.890 

27 

7.354 

23.93 

7979.5 

4371.2 

3.890 

28 

7.782 

24.39 

8033.2 

4342.6 

3.439 

29 

8.292 

24.93 

8097.3 

4319.6 

3.439 

30 

9.223 

25.92 

8199.6 

4271.2 

3.439 

31 

10.182 

26.94 

8298.8 

4222.5 

3.056 

32 

11.148 

27.95 

8358.8 

4182.8 

2.732 

33 

12.088 

28.94 

8417.3 

4146.7 

2.732 

34 

14.966 

32.00 

8651.5 

4028.0 

2.297 

35 

17.875 

34.95 

8942.1 

3919.1 

1.965 

36 

20.789 

38.00 

9017.3 

3821.3 

1.711 

37 

23.756 

41.00 

9123.0 

3751.5 

1.550 

38 

26.706 

44.00 

9226.7 

3698.2 

1.471 

39 

29.675 

47.00 

9307.7 

3645.8 

1.370 

40 

31.005 

48.33 

9346.7 

3621.9 

1.347 
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VELOCITY (ft/sec) 




Figure 13. MNASA Nozzle Boundary Layer Edge Pressure 
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TEMPERATURE (*R) 



Figure 14. MNASA Nozzle Boundary Layer Edge Temperature 



Figure 15. MNASA Nozzle Cutoff Analysis Grid 
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Figure 16. Boundary Layer Edge Velocities for the Submerged Region 

Results from two grids were combined in defining the heating conditions 
throughout the nozzle. As stated above, the MEIT module was used to calculate the 
convective heating conditions. Results from these calculations are presented in Figure 17. 
This figure presents the convective transfer coefficient as a function of streamlength. It 
appears in this figure that the results from the cutoff grid are shifted slightly to the left. 
This is due to a slight difference in the streamlength which is calculated for the two grids. 
Figure 18 presents results as a function of axial position in the nozzle. The zero axial 
position is taken as the nozzle throat. The results are presented beginning at the nose of 
the nozzle (i.e., the point located furthest into the motor). The shift is no longer evident in 
this figure. Results aft of the throat agree precisely for the two grids. However, the 
cutoff grid calculates slightly lower heating rates in the throat region. Also, from Figure 
17 , the cutoff grid yielded higher heating rates in the submerged region of the nozzle, 

(i.e., streamlengths less than 9 inches). 
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HEAT TRANSFER COEFFICIENT (LBM/ft*-sec) HEAT TRANSFER COEFFICIENT (LBM/fH-sac) 



Figure 17. MNASA Nozzle Convective Heating Conditions 



Figure 18. MNASA Nozzle Convective Heating Conditions as a Function of Axial 

Position 
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The sensitivity of the convective heating conditions, in the submerged region of the 
nozzle, to the fluctuating flow were also investigated. Results for this comparison are 
shown in Figure 19. This figure shows the convective transfer rate as a function of cycle. 
These results are taken from the calculation with the cutoff grid. The difference in the 
calculated rates is very small. Therefore it was decided to use the results from the 3000 
cycle solution as the steady state solution. A summary of the nozzle heating conditions is 
contained in Table 4. This table presents the radiant heat flux, Q rac j, the radiation view 
factor, F, the recovery enthalpy, H r , the convective transfer coefficient h/Cp, and the edge 
pressure, P e , as a function of nozzle station. The location of each of these stations is 
shown in Figure 20. Nozzle erosion measurements have been performed at each of these 
stations. Erosion calculations (see Section 4) were also performed at these same 
positions. The radiant heat flux and the radiation view factors given in the table were 
determined using empirical techniques available in the NAT code. It should be mentioned 
that this table also includes interface heating effects which are described in the upcoming 
section. 

The particle radiation is calculated with NAT assuming a parallel plate model. 

This model applies to aluminized propellants and assumes that the particle -laden stream of 
combustion products is optically thick and exchanges radiant energy with the surface as if 
the stream and wall were parallel plates. In this way, multiple reflections between wall and 
stream are considered. In addition, both the stream and wall were assumed to behave as 
gray bodies which emit and reflect radiant energy diffusely. Based on the above 
assumptions, the net radiant heat flux is given as: 

q n «nd = ee ffa (T^-T 4 w) (1) 


where. 


£eff 

8s 

£w 

a 

Ts 

Tw 

C 

n 

P 

D 


effective emmissivity = 

Vew+Ve.-1 


particle-laden stream emissivity = l-exp( -C-^-pD 

V 16 


wall emissivity 

Stefan-Boltzmann constant 

freestream (edge) temperature 

wall temperature 

empirical constant = 0.808 

percentage of aluminum loading 

local density of propellant combustion species (lb/ft 3 ) 

local beam length, usually taken as the diameter (inch) 


Table 3 presents quantities Qrad an ^ 8 which are input into the CMA thermal 
analysis code. These quantities are given by: 
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HEAT TRANSFER COEFFICIENT (LBM/fP-sec) 


Qrad=^— 

£ w 

F = — 

£w 

The convective transfer coefficients presented in the figures and table were 
determined based on the average chamber pressure for the 48-5 firing. Convective 
transfer rates for a given point in time, t , were determined with the following expression 


(h/Cp) t J P t V : 
(h / Cp) avg y Pavg J 



Figure 19. Submerged Region Convective Heating Conditions, Cutoff Grid 
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Figure 20. MNASA Nozzle Thermal Analysis Locations 
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Table 4: MNASA Nozzle Heating Conditions 


■■■■ 

Qrah - 

MmxitffimM 

j p 


§J fc/Cp III# 

(Ibm/ft^-sec) 

i Pe(atm) 

4 

798.92 

1.000 

1160.00 

0.041 

44.900 

5 

798.92 

1.000 

1160.00 

0.065 

45.000 

6 

797.57 

1.000 

1160.00 

0.069 

44.700 

7 

796.27 

1.000 

1160.00 

0.078 

44.600 

8 

794.08 

1.000 

1150.00 

0.087 

44.400 

9 

800.31 

1.000 

1160.00 

0.136 

45.000 

10 

803.52 

1.000 

1170.00 

0.186 

45.300 

11 

806.68 

1.000 

1180.00 

0.299 

45.700 

12 

795.13 

1.000 

1180.00 

0.436 

45.700 

13 

730.62 

0.960 

1 140.00 

0.707 

40.900 

14 

677.90 

0.935 

1120.00 

0.883 

37.500 

15 

657.60 

0.920 

1120.00 

0.903 

36.600 

16 

555.44 

0.870 

1 1 10.00 

1.050 

29.700 

17 

468.98 

0.835 

1065.60 

0.975 

26.280 

18 

433.77 

0.831 

1070.00 

0.936 

19.800 

19 

293.71 

0.754 

1020.00 

0.695 

12.060 

20 

190.12 

0.657 

878.80 

0.416 

7.281 

21 

157.57 

0.603 

886.76 

0.381 

6.208 

22 

135.90 

0.563 

882.14 

0.353 

5.507 

23 

116.89 

0.521 

878.96 

0.328 

4.923 

24 

101.98 

0.486 

874.54 

0.309 

4.458 

25 

93.80 

0.466 

872.51 

0.297 

4.191 

26 

84.79 

0.444 

870.51 

0.500 

3.890 

27 

84.79 

0.444 

870.51 

0.460 

3.890 

28 

72.07 

0.412 

865.24 

0.417 

3.439 

29 

72.07 

0.412 

865.24 

0.409 

3.439 

30 

72.07 

0.412 

865.24 

0.372 

3.439 

31 

61.71 

0.384 

860.50 

0.343 

3.056 

32 

53.29 

0.366 

845.45 

0.319 

2.732 

33 

53.29 

0.366 

845.45 

0.301 

2.732 

34 

43.25 

0.332 

841.47 

0.266 

2.297 

35 

36.05 

0.310 

839.64 

0.228 

1.965 

36 

30.56 

0.290 

836.46 

0.199 

1.711 

37 

27.42 

0.275 

830.37 

0.180 

1.550 

38 

25.98 

0.265 

825.20 

0.167 

1.471 

39 

24.65 

0.253 

810.23 

0.151 

1.370 

40 

24.38 

0.250 

804.47 

0.145 

1.347 


asrm 
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3.1.1 Material Interface Heating 

It is commonly known that relatively small surface discontinuities can modify 
flowfields and augment convective transfer rates. A surface discontinuity is known to 
occur in the nozzle exit cone at the juncture of the standard and low density carbon 
phenolic materials. This surface discontinuity is the result of differential ablation rates for 
the two materials. The low density material ablates at a faster rate than the standard 
density material and thus a rearward facing step is formed. This step results in separated 
flow and boundary layer reattachment just aft of the material interface. Initial attempts 
were made to model this flow with the KIVA module of NAT in order to determine 
augmented heating rates. These calculations used the post test nozzle shape, i.e., eroded 
surface contour, but never yielded satisfactory results. Therefore, it was decided to 
address this issue using simplified procedures. 

A literature search was performed and yielded a number of applicable references, 
see References 3 through 6 . After reviewing these documents it was decided to simply 
calculate the heating conditions in this region by assuming that the boundary layer 
reattached just downstream of the material interface. This calculation was performed by 
using the flow results for the uneroded nozzle contour (described in Section 3.1) and 
assigning the boundary layer starting point as the material interface location. Results from 
this calculation are shown in Figure 21. This figure compares the convective transfer 
coefficients for the unperturbed and augmented heating conditions. The streamlength in 
this case is measured from the material interface location. The figure presents results 
beginning at this location and continuing to the exit plane of the nozzle. Starting the 
boundary layer at the interface position drastically increases the convective transfer rates, 
(approximately 15% nearest the interface). However, this augmentation rapidly falls off 
and there is only a slight augmentation effect at the exit plane location. How these results 
are used in subsequent thermochemical erosion calculations will be further described in 
Section 4. 


3.1.2 Particle Flow Calculations 

Two phase flow calculations were performed for the MNASA motor with the 
KIVA module of NAT. These calculations assumed that all of the aluminum in the 
propellant reacted to form liquid AI 2 O 3 . Several initial calculations were performed for 
the nozzle considering 8 and 10 qm diameter particles. These calculations were found to 
be very time consuming and additionally no particle impacts were calculated to occur with 
the nozzle surface. Based on these initial experiences a less refined grid was developed to 
decrease computation time. This grid is shown in Figure 22. To perform these 
computations particles are modeled in groups called parcels. These parcels are 
continuously injected into the flow as indicated in the figure. Prior to injecting the 
particles, a steady state 2-D gas flow solution was obtained. First obtaining this steady 
state gas flow solution is not explicitly required for the particle flow calculation, but was 
simply the method selected for this study. This method was selected because the KIVA 
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module is a time accurate Navier Stokes solver. Because of this, it was desired not to 
include the transient effects in going from the initial flow conditions to the converged gas 
flow conditions. 
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h/Cp (Ibm/ft2-sec) 



Figure 21. Augmented Heating Rates at the Exit Cone Material Interface 



Figure 22. MNASA Particle Flow Grid 


2 ‘ 



Based on data of Salita (Reference 7) a calculation was performed for a 120 pm 
particle diameter. Salita performed experiments which showed that in a motor chamber a 
significant amount of particles with diameters greater than 100 pm are formed. Particle 
size data from Salita is shown in Figure 23. This figure shows the mass fraction of 
particles smaller than a given size. Both the data and Salita’s curve fit are presented. 
Notice that about 65% of the particle mass consists of particles with diameters less than 
10 pm. Nearly all of the remaining particle mass , i.e., 35%, is contained in particles with 
diameters greater than 100 pm. This data is taken from an experiment with a 500 psi 
chamber pressure. The significant quantity of large particle diameters may be inconsistent 
with particle sizes which are traditionally considered in solid rocket motor nozzle analysis. 
However, this data is consistent with other studies (References 8) which show that large 
particles can exist in the motor chamber but will breakup in the exit cone and exit the 
nozzle much smaller in size. 



Figure 23. Particle Size Distributions from Reference 7 


25 



Results from the 120 pm particle flow computation are presented in Figures 24 
through 28. The total particulate mass flow rate for the motor is determined from one 
dimensional isentropic flow solutions. The total number of parcels used for the 
computation is then selected. After determining the mass of each particle it is then 
possible to determine the number of particles per parcel and the required injection rate. 

For this calculation each parcel represents 18,729 particles and parcels are injected into 
the flow at a rate of 748,000 parcels/second (this rate is for a 0.5 degree sector of the 
motor). The injection was performed in such a fashion that a uniform particle mass flux 
was entering the flow domain. Because of this, a higher parcel injection rate is required at 
the larger radial positions (because of the larger area in these regions). These figures 
present parcel velocity plots as a function of time. Results are presented for only selected 
parcel numbers; 1-100, 400-500, 900-1000, 1900-2000, and 2900-3000. Following are 
comments relative to these results: 

1. In the entrance region of the nozzle, particle velocities are greater at locations away 
from the nozzle centerline. 

2. Particle impacts are calculated to occur for both the aft closure and submerged 
portions of the nozzle. 

3. Particles are accelerated radially inward in the converging portion of the nozzle and 
tend to remain very tightly grouped near the nozzle centerline through the exit cone. 

Figure 29 presents parcel velocities plotted for every fifth parcel. This plot helps to show 
the trajectory of each parcel since they tend to follow one after another. Also identified in 
this figure is the region within the nozzle for which particle impacts are calculated. Since 
parcels are injected at discrete locations, impacts also occur at discrete locations. 

However, in reality these impacts cover an area just as the particle inflow, although 
modeled with discrete injection points, covers as area. The impact conditions at location 1 
(identified in the figure) are a velocity of 140 ft/sec and an impact angle of 46 degrees. 

The conditions at location 2 are 154 ft/sec and an impact angle of 71 degrees. 

The mass flux of the particles was determined for the two locations by summing 
the total number of parcel impacts over time and knowing the mass of each parcel. 
Additionally, the area of impact was taken as shown in the figure. With this approach, the 
total particle mass fluxes were found to be 12.4 (lbm/ft^-sec) and 15.7 (Ibm/ft^-sec) at 
locations 1 and 2, respectively. Of course, this assumes that all particles are 120 |J.m as 
used in this analysis. The data from Salita"7 suggests that the 120 pm diameter particles is 
probably representative of 35% of the particle mass. Based on this, it was assumed that 
only 35% of the mass flux actually was impacting the nozzle (previous calculations for 10 
pm particles had shown no impacts). With this, the mass fluxes become 4.34 (lbm/ft^-sec) 
and 5.495 (lbm/ft--sec) for the two locations. This method of applying particle size 
effects was used in lieu of attempting flow computations with particle size distributions. 
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The next task was the determination of impact angle and velocity for each of the 
nozzle measurement positions for which particle impact occurred. The curves in Figure 
30 were used to determine this information. This figure presents particle velocity, particle 
direction and body angle as a function of axial position in the nozzle. The body angle is 
the angle between the local nozzle surface tangent and the nozzle centerline. The body 
angle plot is simply determined from the nozzle geometry. The velocity plot was 
constructed by using the velocities at the two impact locations and linearly extrapolating. 
The particle direction plot was constructed by using the two impacting trajectories and 
inferring the trajectories for the outer regions of the impacts. These trajectories were 
inferred based on both the impacting and non impacting particle trajectories. With this 
data, it was then possible to determine the particle impact conditions for each of the 
analysis positions. This information is summarized in Table 5 for measurement locations 
where particle impacts are calculated to occur. The impact information contained in Table 
4 was used in subsequent material erosion calculations which are described in Section 4. 

It should be noted that the flexseal specimen used in these tests is located at axial positions 
between -0.7 inch and 3.55 inches. Obviously, this specimen is experiencing particle 
impact as well as convective and radiative heating. 
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Figure 30. Particle Impact Conditions for the MNAS A Nozzle 
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Table 5. MNASA Nozzle Particle Impact Conditions 


Station 

Axial l! 
Position 

lliin.) 

§f Radial: ^ | 
Position 

|Surfiacti| 

:\Ahgle|f 

: (deg) 

Diretfioii 
- j ‘ tfeg) 

Impact 
i; Angle 
(deg) 

Velocity 
; (ft/sec) 

!M . nr :t!|I 

(lbm/ft 2 

|@ sec) ' ' 

6 

1.393 

10.135 

14 

m 

35 

126 

4.34 

7 

0.365 

9.863 

14 

E9 

38 

133 

4.34 

8 

-0.701 

9.590 

14 

149 

47 

139 

4.34 

9 

-1.923 

9.307 

14 

140 

54 

147 

4.34 

10 

-2.863 

9.044 

17 

131 

66 

152 

5.495 

11 

-3.815 

8.59 

38 

118 

100 

157 

5.495 

12 

-4.232 

7.694 

80 

111 

149 

161 

5.495 


3.2 ASRM NOZZLE FLOWFIELD AND HEATING RATE CALCULATIONS 

Nozzle flowfield, heating rate , and surface thermochemistry calculations have 
been performed for the ASRM nozzle using the NAT computer code and the same 
procedures employed for the MNASA nozzle. This motor uses the same propellant as 
utilized in the MNASA motor and therefore the propellant gas formulation given in Table 
2 was again used. Figures 3 1 and 32 present two different analysis grids used to calculate 
the flow within the ASRM nozzle. Figure 31 is the analysis grid for time zero, i.e., with 
the propellant grain in its initial configuration. This figure also identifies the propellant 
grain boundary. The inflow conditions were again applied at the left hand side of the 
analysis grid. Figure 32 presents the analysis grid used for the nozzle at the final point in 
time and will be referred to as the end of bum (EOB) grain. It should be pointed out that 
the motor actually attains this configuration at about 112 seconds into the bum. The aft 
motor section is burned out at this point but the forward section has propellant remaining 
and continues to bum. The gridding of these two configurations is similar to the MNASA 
motor in that grid refinement is greatest in regions where the largest gradients are 
expected. The average chamber pressure of 616 psi was used for both geometries. 

Because of the differences in the motor grain geometries, the inflow velocities for 
the two configurations are drastically different. The initial grain configuration has an 
inflow velocity of about 1050 ft/sec whereas the EOB grain has an inflow velocity of only 
290 ft./sec. It was expected that this difference in inflow conditions could potentially 
cause significant differences in the nozzle flowfields and heating conditions, especially in 
the submerged and the throat entrance regions of the nozzle. 
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Flow results for the initial grain geometry are presented in Figures 33 through 36. 
These are the steady state flow results obtained after a 2800 cycle calculation. Figure 33 
presents flow velocity vectors. Figures 34 through 36 present pressure, temperature, and 
Mach contours. The most significant feature of these results is the very large gradients 
which exist just downstream of the nozzle throat. 

Flow results for the EOB grain geometry are presented in Figures 37 through 40. 
These results required 2000 cycles to reach steady state conditions. Besides the obvious 
differences in the entrance and submerged regions of the nozzle, there are only very small 
differences in the flow results compared to the initial grain geometry. A precise 
comparison is not easily made because of the slight differences in the value of the contour 
lines. However, the features and shape of the contours are very similar. This is most 
evident in the forward exit cone region. The slight differences in the maximum and 
minimum temperatures and pressures are most likely just a result of the degree of 
convergence in the solution. 

Boundary layer solutions were performed for the nozzle using the same procedures 
described in Section 3.1. Figure 41 shows the location within the nozzle chosen as the 
boundary layer starting point. This location was selected based on the flow direction 
calculated for this area, i.e., this appeared to be the initial point where the flow followed 
the nozzle surface continuously until exiting the nozzle. The heating rates in the 
converging/diverging portions of the nozzle were found to be insensitive to this starting 
point. 



Figure 33. ASRM Nozzle Velocity Vectors, Initial Grain 
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Figure 34. ASRM Nozzle Pressure Contours, Initial Grain 



Figure 35. ASRM Nozzle Temperature Contours, Initial Grain 
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Figure 38. ASRM Nozzle Pressure Contours, EOB Grain 



Figure 39. ASRM Nozzle Temperature Contours, EOB Grain 
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M max = 2.74 
Mmin = 0.152 
AM = 0.152 

Figure 40. ASRM Nozzle Mach Contours, EOB Grain 



Figure 41. Boundary Layer Starting Point for ASRM Nozzle 
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The boundary layer edge conditions and heat transfer coefficients for the nozzle 
are presented in Figures 42 through 49. Results are presented in two groups. The first 
group. Figures 42 through 45, covers positions within the nozzle beginning at the nose of 
the nozzle and continuing downstream through the exit cone. The second group. Figures 
46 through 49 , begins at the nose and runs into the submerged region of the nozzle. 

These figures present velocity, pressure, temperature and convective heat transfer 
coefficient as a function of axial position within the nozzle. These axial positions are 
taken from the drawing for the nozzle. The zero point is located at the head end of the 
motor. The nose of the nozzle is located at an axial position of 421.8 inches. The nozzle 
throat is located at 446.6 inches. Following are comments relative to these results: 

1. There is a sudden expansion of the flow just aft of the nozzle throat. This can be 
seen in Figures 42 through 44 but is most pronounced in Figure 45, the heat transfer 
coefficient. In this figure, the heat transfer coefficient reaches a maximum just 
upstream of the throat and then drops off very suddenly. 

2. Temperatures and pressures in the submerged region are nearly constant. 

3. Flow velocities in the submerged region of the nozzle are about 100-150 ft/sec in the 
aft regions (nearest the nose) but drop off to about 50 ft/sec in the forward region. 
The heat transfer coefficient follows the same trend. It has a value of about 0. 1 
(lbm/ft^-sec) and then drops off to about 0.025 (lbm/ft^-sec). 



Figure 42. ASRM Nozzle Edge Velocity 
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Figure 43. ASRM Nozzle Edge Pressure 



Figure 44. ASRM Nozzle Edge Temperature 
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Figure 46. ASRM Nozzle Edge Velocity, Submerged Region 
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Figure 47. ASRM Nozzle Edge Pressure, Submerged Region 
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Figure 48. ASRM Nozzle Edge Temperature, Submerged Region 
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420.0 425.0 430.0 435.0 440.0 445.0 450.0 

AXIAL POSITION (in) 

Figure 49. ASRM Nozzle Edge Heat Transfer Coefficients, Submerged Region 

Results which have been presented thus far are for the initial grain configuration. 

As previously mentioned, the flow results for the EOB configuration were relatively 
similar. Just how similar can be seen in Figures 50 and 51. Figure 50 shows the heat 
transfer coefficients calculated fer ]%%~%%or the nozzle downstream of the nose. This figure 
compares the results for the initial and end of bum grain geometries. The initial grain 
geometry appears to have slightly higher convective heating rates in the nozzle entrance 
region but otherwise the heating rates are nearly identical. Similar results are presented in 
Figure 51 for the submerged portion of the nozzle. Not too surprisingly, the heating 
conditions are somewhat different. The EOB configuration tends to experience lower 
heating conditions near the nose of the nozzle (i.e., axial positions around 425). 

Several locations within the ASRM nozzle were selected for subsequent thermal 
and erosion calculations. These locations are identified in Figures 52 and 53. Figure 52 
shows the analysis locations in the entrance region of the nozzle and Figure 53 identifies 
the positions in the exit cone region. Table 6 presents the heating conditions calculated 
for each of these analysis positions. All of the heating conditions presented are for the 
initial grain configuration. Only for the submerged region, station 1, were time effects 
included in the erosion calculations. This table also includes interface heating conditions 
(beginning at Station 11) which will be discussed in the following section. 

One area of the nozzle receiving special attention was the flex seal. Tests 
conducted with the MNASA motor had shown the flex seal specimen to experience a 
significant amount of erosion. Therefore, accurately defining the heating conditions for 
the ASRM flex seal was of increased importance. Table 7 presents a summary of the 
results calculated for the flex seal region. Somewhat surprisingly, the conditions are 
expected not to change very much during the bum. Gas velocities adjacent to the flex seal 
range from about 30 ft/sec to slightly over 100 ft/sec. The convective heating rate on the 
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flex seal is about 0.06 lbm/ft 2 -sec. This is only about 25% less than the conditions 
calculated for the MNASA motor. 



h/Cp (Ibm/ft2-sec) h/Cp (Ibm/ft2-sec) 









Figure 52. ASRM Nozzle Thermal Analysis Locations, Forward Region 




Figure 53. ASRM Nozzle Thermal Analysis Locations, Aft Region 



Table 6. ASRM Nozzle Heating Conditions 


Hi iv 


mmm 



752.81 


792.64 


794.68 


751.00 


700.60 


653.94 


519.34 


322.21 


282.70 


213.24 


176.08 


127.22 


87.36 


83.73 


72.83 


71.22 


68.81 


0.998 


0.989 


0.964 


0.942 


0.898 


0.841 


0.833 


0.792 


0.783 


0.769 




1122.2 


1186.4 


1199.6 


1191.2 


1 177.4 


1172.5 


1175.7 


882.3 


865.6 


854.0 


851.0 


842.9 


803.9 


800.4 


800.5 


800.5 


798.3 



0.335 


0.260 


0.370 


0.221 


0.145 


0.122 


0.116 


0.109 
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Table 7. AS1 

TM Flex Seal Heating Conditions 

Pressure (atm) 

iBMI 

H/Cp ; 

iiii 


2 -sec): 

40.7 

1188 

0.0687 (IG) 
0.0608 (EOB) 

1.0 

752 


The ASRM Motor experiences a relatively large variation in pressure during the 
course of the motor bum, see Figure 6. Because of this, the approximation given in 
equation 1 to account for chamber pressure variations was compared to accurate 
boundary layer solutions for another pressure condition. Table 8 presents results from this 
comparison. This table presents the heat transfer coefficient calculated by NAT for two 
different chamber pressures, 616 psi and 910 psi. Both of these calculations were 
performed for the initial grain geometry. Results are presented at five locations within the 
nozzle. The pressure of 910 psi was selected because it is the maximum pressure the 
ASRM motor will experience. The heat transfer coefficient calculated with equation 1 is 
shown for comparison and the percent difference is also given. Equation 1 does an 
excellent job of predicting heating conditions near the nozzle throat but then slightly 
underpredicts heating in the exit cone. The differences, however, were considered small 
enough that equation 1 was used for locations throughout the nozzle. 


Table 8. Comparison of Calculated and Approximated Heat Transfer Coefficients 


Location ; 

Axial Position 
tins’) 

. (hM616 

flhni/ft2-s€C> 

m 

<PRAT) ,S (h/cp)616 

llbm/rt 2 -secl 

difference 

Throat 

Entrance 

425.1 

0.5820 

0.7928 

0.7948 

.25 

Throat 

442.0 

0.7818 

1.067 

1.0677 

.07 

Entrance 





-1.61 

Throat 

446.6 

0.6590 

0.9147 

0.9000 


FVVD Exit 
Cone 

454.9 

0.3543 

0.4993 

0.4839 

-3.08 

AFT Exit 
Cone 

489.2 

0.2194 

0.3043 

0.2996 

-1.54 
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3.2.1 Material Interface Heating 

A material interface heating calculation was performed for the ASRM nozzle using 
the same procedures as described in Section 3.1.1 for the MNASA nozzle. This 
calculation was performed for the location in the nozzle where the material changes from 
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standard density MX4996 carbon phenolic to low density FM5939 carbon phenolic. The 
"interface heating" condition was determined by simply using the flow results for the initial 
grain geometry and starting the boundary layer solution at the material interface location. 
Results form this calculation are presented in Figure 54. This figure compares the 
augmented heating results (i.e., starting the boundary layer solution at the material 
interface) to the unperturbed solution. The results are presented as a function of 
streamlength beginning at the interface and continuing to the exit plane of the nozzle. 
Heating conditions nearest the interface are increased by nearly 100%. However, this 
effect falls off down the exit cone until the augmentation is only about 10% at the exit 
plane of the nozzle. 



3.2.2. Particle Flow Calculations 

Particle flow calculations were performed for the ASRM nozzle using the same 
procedures as for the MNASA nozzle. These calculations incorporated lessons learned 
from the MNASA calculations to improve upon accuracy and efficiency. The calculations 
were performed to determine the location and severity of particles impacting on the nozzle 
surface. To perform these calculations, three particle flow grids were generated. These 
grids, shown in Figures 55 through 57, are somewhat more coarse than the grids used for 
the flow solutions. Some degree of grid refinement was maintained in the regions of the 
nozzle adjacent to the nose. This was done so that particle impacts in this region could be 
accurately calculated. The three grids represent three separate points during the motor 
burn. The first grid is for the initial grain configuration. The second grid is for an 
intermediate point in time, 24 seconds into the firing. The third grid is the EOB 
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configuration and represents the motor geometry 112 seconds into the firing. Four 
different particle diameters (1, 3.5, 10, and 120 Jim) were considered in these calculations. 
Table 9 presents a summary of the various particle flow calculations which were 
performed. 
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Figure 56. ASRM Nozzle Particle Flow Grid, Intermediate Grain Configuration 
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Figure 57. ASRM Nozzle Particle Flow Grid, EOB Grain Configuration 


Table 9. ASRM Nozzle Particle Analysis Cases 
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Figures 58 through 66 present parcel velocity results for the initial grain 
configuration. Results are presented for particle sizes of 1,10 and 120 pm. Three figures 
are presented for each particle size. Each figure represents a slice in time and shows 
velocities for selected parcels. The sequence of figures shows how the particle flow 
develops. Notice that it takes less than 0.02 second for the particle flow fields to fully 
develop. Figures 67 through 69 present similar results for all three particle sizes with the 
EOB configuration. For this configuration, about 0.04 second is required for a fully 
developed particle flowfteld. This is primarily because of the much lower inflow velocities 
for this configuration. These figures present results for a single point in time and show the 
fully developed particle flow. Following are comments relative to the figures: 

1 . The 1 and 10 pm panicles tend to follow the flow stream lines. They both come 
very close to impacting the nozzle in the throat entrance region but turn prior to 
impacting. 

2. The 120 pm particles do not follow the flow stream lines as closely and impact the 
nozzle surface in the throat entrance region early in time and in the submerged region 
late in time (i.e., at end of bum). 

3. The smaller particles (i.e., the 1 and 10 pm diameters) follow the flow stream lines 
and expand to near the surface in the nozzle exit cone. There is very little difference 
between the initial and EOB configurations 

4. The 120 pm particles tend to form more of a core in the exit cone and do not even 
come close to the exit cone surface. This core becomes tighter later in time, i.e., is 
located nearer the centerline. 

Based on the diversity of the impact results for the initial and EOB configurations, 
it was determined that it would be necessary to run particle flow calculations for at least 
one additional point in time. This was necessary so that particle impact conditions could 
be defined as a function of time for the erosion calculations. The grain configuration 
shown in Figure 56, corresponding to 24 seconds after motor ignition, was selected as the 
intermediate point. Calculations were performed for only the 120 pm particle size. 

Results for this calculation are presented in Figure 70. It appears that this calculation does 
a good job of bridging the impact conditions form the initial and EOB grain 
configurations. Impacts are calculated to occur in both the throat entrance and submerged 
regions of the nozzle. 

One final particle flow calculation was performed for the nozzle. This calculation 
was performed explicitly to try to produce impacts in the exit cone. This calculation was 
performed for the initial grain geometry and a 3.5 pm particle diameter. This particle 
diameter was selected based on information from Eric Wemimont of Aerojet^- Results 
from this calculation are given in Figure 71. Particle trajectories are calculated to come 
slightly closer to the exit cone than for either the 1 or 10 pm particles, but still no impacts 
are calculated to occur. 
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Figure 58. ASRM Particle Flow Results, Initial Grain 1 Jim Particles. 

Time = 0.0090 second 



Figure 59. ASRM Particle Row Results, Initial Grain 1 |im Particles, 

Time = 0.01 16 second 



Figure 60. ASRM Particle Row Results, Initial Grain 1 |im Particles, 

Time = 0.0343 second 


Figure 61. ASRM Particle Flow Results, Initial Grain 10 Particles, 

Time = 0.0086 second 



Figure 62. ASRM Particle Flow Results, Initial Grain 10 |im Particles, 

Time = 0.0104 second 



Figure 63. ASRM Particle Flow Results, Initial Grain 10 |J.m Particles, 

Time = 0.0195 second 



Figure 64. ASRM Particle Row Results, Initial Grain 120 |im Particles 

Time = 0.0082 second 



Figure 65. ASRM Particle Flow Results, Initial Grain 120 |im Particles, 

Time = 0.0110 second 



Figure 66. ASRM Particle Row Results, Initial Grain 120 |im Particles, 

Time = 0.0163 second 
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Figure 67. ASRM Particle Flow Results, End of Bum 1 Jim Particles 



Figure 68. ASRM Particle Flow Results, End of Bum 10 ^im Particles 



Figure 69. ASRM Particle Flow Results, End of Bum 120 (im Particles 




Figure 70. ASRM Particle Flow Results, Intermediate Grain 120 (am Particles 



Figure 71. ASRM Particle Flow Results, Initial Grain, 3.5 Jim Particles 

Particle mass fluxes were determined for the ASRM nozzle using the same 
procedures as for the MNASA nozzle. The total number of parcel impacts was 
determined for a period of time and then the particle mass flux per unit area was 
determined. The Salita^ data was again used to account for particle size distributions. 
Based on Figure 23 the impacting mass flux was taken as 35% of the calculated amount 
for the 120 [im particles.. These mass flux results are presented in Figure 72. This figure 
shows the calculated mass fluxes as a function of radial position in the nozzle. Results are 
presented for three points in time, initial (INTT), intermediate (INT), and end of burn 
(EOB). For reference, the throat radius is 27.2 inches and the nose radius is 40.07 inches. 
The corresponding impact angle and impact velocity results are presented in Figures 73 
and 74. Following are comments relative to these figures. 

1. The greatest mass fluxes occur early in time and are approximately five times greater 
than those generated in the MNASA motor. 

2. The impacting mass fluxes drop off significantly with time and for the EOB 
configuration are actually less than the MNASA levels. 
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3. For the initial grain, the impact angles are relatively low near the nozzle throat and 
approach 70 degrees near the nose. The same is true for the intermediate grain. 

4. Impact angles for the submerged region are very small for the intermediate time and 
approach 90 degrees for the EOB configuration. 

5. Impact velocities are highest early in time and drop off to around 100 ft/sec at the 
end of bum. The higher velocities are a result of the high inflow velocities which are 
caused by the relatively small difference between the throat area and grain bore area 
early in time. These inflow velocities drop off significandy as the grain recedes. 
Impact velocities drop off a corresponding amount. 



Radius (in) 


Figure 72. ASRM Particle Mass Fluxes 

The data in Figures 72 through 74 was used to determine particle impact 
conditions at each of the affected ASRM nozzle analysis locations. The calculations 
showed impacts to occur only at analysis locations 1 through 5. These conditions are 
presented in Figures 75 through 77. These figures present impacting mass flux, incidence 
angle and impact velocity as a function of time. 
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Impact Angle (deg) Mass Flux (Ibm/ft2-sec) 
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Figure 75. Particle Mass Fluxes at ASRM Analysis Locations 



Figure 76. Particle Impact Angles at ASRM Analysis Locations 
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Figure 77. Particle Impact Velocities at ASRM Analysis Locations 


3.3 SURFACE THERMOCHEMISTRY 

Surface thermochemistry for each of the nozzles was calculated with the ACE 
module of NAT. These calculations were performed based on the assumption of chemical 
equilibrium between the pyrolysis gases, the charred surface material and the propellant 
combustion gases. The composition of the propellant combustion gases was based on an 
isentropic flow chemical equilibrium composition at the location of interest. The 
composition of the pyrolysis gases for each of the materials is given in Appendix A. 
Figures 78 through 80 present representative results for each of the materials considered. 
These results are presented only for the MNASA motor since the propellant compositions 
of the two motors are identical. Each figure presents nondimensional surface removal 
rates, B' c , as a function of surface temperature. Results are presented for a range of 
nondimensional pyrolysis gas rates, B'g All of these results are typical of carbon 
phenolics. There is a diffusion limited plateau up to temperatures of about 5000 R at 
which point the B' c values begin to increase rapidly with temperature. Larger values of 
B'g tend to reduce the level of the plateau and therefore result in lower B' c 's This 
information along with corresponding surface energy fluxes is supplied to the thermal 
analysis codes and used in their surface energy balance calculations to determine 
thermochemical erosion rates. 
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Figure 78. Non Dimensional Ablation Rates for FM5055 Carbon Phenolic 



Figure 79. Non Dimensional Ablation Rates for MX4996 Carbon Phenolic 
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4.0 NOZZLE MATERIAL RESPONSE CALCULATIONS 


This section presents results from material response calculations which have been 
performed for both the MNASA and ASRM nozzles. These calculations include surface 
erosion and in-depth heating of the nozzle materials. The procedures used in these 
analyses have been to predict erosion and char depths for the MNASA nozzle and 
compare to the measured data. This comparison was made to validate the analysis 
procedures and material models which had been used. Calculations were then performed 
for the ASRM nozzle to evaluate the adequacy of the design. 

4.1 ANALYSIS PROCEDURES 

Thermal and ablation calculations have been performed with the Charring Material 
Ablation computer code (CMA, reference 11) and the Aerotherm Axi-Symmetric 
Transient Heating and Material Ablation computer code (ASTHMA, reference 12). CMA 
was used in all nozzle positions where a one dimensional heat conduction response was 
considered a valid assumption. The ASTHMA code was used only to analyze the flex seal 
specimen in the MNASA motor tests (This analysis was performed in conjunction with the 
development of a polyisoprene model and is described in Appendix B). Surface boundary 
conditions and thermochemistry tables were generated with the NAT computer code and 
described in the previous section. The material models which were used in these analyses 
are described in Appendix A. It should be pointed out that the models used for the 
MX4996 and FM5939 carbon phenolics should be considered preliminary since they are 
based on laboratory data and inferred from other carbon phenolic materials. Additional 
testing in the Southern Research Analog Facility and the Aerotherm Arc Plasma Generator 
are needed to fully characterize the material and develop a model which can be used with a 
high level of confidence. This is especially true for the in-depth heating and 
decomposition calculations. The calculated surface erosion rates would not be expected 
to change significantiy with a modified model. 

To account for the augmented heating which occurs at the material interface in the 
exit cone the augmented heating rates presented in Figures 21 and 54 were applied when 
the differential recession was calculated to exceed the local boundary layer thickness. The 
differential recession was calculated based on unaugmented rates. For the MNASA 
nozzle, this condition (i.e., differential recession greater than boundary layer thickness) 
was calculated to occur at 8.8 seconds into the motor firing. The boundary layer thickness 
at the material interface location was calculated to be 0.010 inch. For the ASRM nozzle 
the augmented heating conditions were not applied until 48 seconds into the firing. This is 
because the boundary layer thickness for this nozzle is calculated to be 0.068 inch and 
therefore a much greater amount of time must elapse for the differential erosion to exceed 
the boundary layer thickness. 

A model to predict the surface erosion of the polyisoprene was developed under 
this program. This model was developed to match the surface erosion rates of the 
polyisoprene sample in the MNASA motor. This model treats the material as a surface 
sublimer. A more detailed description of the model is contained in Appendix B. Again, 
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this model is based on a limited amount of data and any calculations for the polyisoprene 
should be considered preliminary. 

An erosion model, described in Reference 9, was used to account for material 
removal due to particle impacts. This model was developed for graphite but was used for 
the carbon phenolic materials in the MNASA and ASRM nozzles. It was necessary to use 
this graphite material model since no particle erosion data exists for carbon phenolic. The 
procedure for modeling the erosion, described in Reference 9, treats material removal due 
to both chemical and mechanical effects of the particulates. Mechanical effects are 
modeled with an expression of the form: 

G me ch = aDbvC(sina)d (3) 


where. 


G 

= ratio of mass of target material removed to 


mass of incident particles 

V 

= particle velocity 

D 

= particle diameter 

a 

= angle of impact relative to the surface 

a,b,c,d 

= empirically determined constants 


The study in Reference 9 developed the following expression for graphite mechanical 
erosion. 

Gmech = 5 .342 x 10 ' 7 D 0 - 799 V 1 - 24 (sin a ) 0 - 299 ( 4 ) 

The units for D and V are pm and ft/sec, respectively. 

The Reference 9 study also found that material removal can occur due to chemical attack 
of either aluminum or alumina (AI 2 O 3 ) particles. This attack was characterized with an 
expression of the form: 


G ch em = a'(DVsina)b' (5) 

where a' and b' are different for aluminum and alumina. The correlation for aluminum 
particles was: 


G c hem = 1.113 x 10' 5 (DVsinoc ) 0 - 8 ( 6 ) 

and for alumina: 

Gchem = 386.7(DVsina)- 1 - 13 ( 7 ) 

These same expressions were used to model the chemical erosion of carbon phenolic, and 
since they are based on relative mass rates (i.e., mass removed/mass incident), they are 
expected to do a reasonable job of predicting erosion for carbon phenolic. This is because 
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the relative mass rates effectively account for density differences between the graphite and 
charred carbon penolic (i.e., this expression will predict greater material removal for 
carbon phenolic than for graphite). For this study, all particles within the two motors 
were assumed to be alumina. This may be realistic for the MNASA motor because of the 
large distance between the propellant and the nozzle ( and therefore sufficient time for 
the aluminum to react with oxygen to form alumina). However, this may not be 
representative for the ASRM where the propellant is packed near the submerged nozzle 
region. 


4.2 MNASA NOZZLE EROSION 

The calculated thermochemical erosion depths for the MNASA nozzle are 
presented in Figure 81. These results are compared to the measured values at locations 4 
through 40. Following are comments relative to these results: 

1. The calculated and measured erosion depths compare very well at locations 4, and 13 
through 33. The exceptions to the good comparison appear to be in regions of the 
nozzle where particle impacts are occurring. 

2. The calculated and measured erosion in the throat region compare exceptionally well. 

3. The calculated erosion levels do a very good job of matching the rapid drop in 
erosion just aft of the nozzle throat (location 18). 

4. The method of accounting for the interface heating appears to do a very good job of 
calculating the increase in erosion at the material interface. It appears to be slightly 
conservative (or this could be partly due to the conductivity model for FM5939, see 
the discussion below) but it does seem to do a good job of matching the 
characteristic of the erosion aft of the interface, i.e., gradually falling off down the 
length of the exit cone. 

Erosion results considering combined thermochemical and particle erosion are 
presented in Figure 82. The particle erosion effects only stations 5 through 12 and 
stations 34 through 40. Therefore, the remaining stations experience only 
thermochemical erosion and the results at those stations are identical to those presented in 
Figure 81. The calculated erosion at stations 5 through 12 compare very well with the 
measured values. Several features of this comparison are important: 

1. The particle impacts are calculated to effect only stations 5 through 12 in the nozzle 
entrance region. This is completely consistent with the measure results. 

2. The magnitude of the erosion compares very well with the measured values. This is 
especially true for stations 5 through 10. At stations 11 and 12 the comparison could 
possibly be much better if the NAT analysis grid had been further refined and the 
impacts at those locations more accurately calculated. The impact conditions were 
extrapolated from the results based on the "area effected method" described in 
section 3. 

3. Based on the relatively good agreement with measured data no attempt was made to 
update or modify the G law correlations from Reference 9. 
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The matching of the particle erosion effects in the nozzle exit cone was based on 
backing out particle impact conditions. Recall from section 3 that no particle impacts 
were calculated to occur in the nozzle exit cone. Therefore, the impacting mass flux was 
backed out to match the erosion data. The impacting particle diameter was assumed to be 
3.5 |im. Since this is such a small particle size the particle velocity was assumed to be 
equal to the local gas velocity. The impacting angle was assumed to be 5 degrees. Based 
on all of these assumptions, the impacting mass fluxes, shown in Figure 83, were 
determined. This figure shows the particle mass fluxes as a function of streamlength. The 
large values of 4-5.5 lbm/ft 2 -sec occur in the nozzle entrance region. Mass fluxes in the 
exit cone are calculated to be only about 0.2 lbm/ft 2 -sec and increase slightly nearest the 
exit plane of the nozzle. 

Figure 84 presents the total affected depth (TAD) calculated for the MNASA 
nozzle and compares it to measured results. TAD is the sum of the erosion depth and the 
post test measured char depth. A summary of the measured and calculated erosion and 
TAD results is also presented in Table 10. TAD’s are presented at the end of the motor 
burn and following two minutes of cool down. This two minute cool down allowed for 
heat to continue to soak into the material and additional charring to occur. The cool 
down environment was not known and therefore it is assumed that only a small amount of 
heat leaves the nozzle surface due to radiation. A radiation view factor of 0.05 was used 
at all stations within the nozzle. 

The TAD predictions for stations 4 through 24 do a good job of bounding the 
data. This is the best that can be expected based on the unknown cool down environment. 
However, beginning at station 26 (this is the location where the FM5939 material begins) 
the predicted TAD's are below the measured values even with the cool down effects. This 
is probably due to the material model for FM5939. The conductivity model for the 
MX4996 was modified (i.e., the conductivities were increased from the laboratory 
measured values) to be similar to other standard density carbon phenolic material models. 
This was not done for the low density FM5939 carbon phenolic. Based on these results 
more data is needed to do a reasonable job of predicting the char depth for the FM5939. 
Based on the TAD results the thermal conductivity would be expected to be greater than 
the values used in the analysis. It is interesting to note that an increase in conductivity 
would result in reduced levels in the calculated thermochemical erosion. This would 
improve the agreement between measured and calculated erosion for the FM5939, see 
Figures 81 and 82. 
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Figure 81. Calculated Thermochemical Erosion Depths for the MNASA Nozzle 



Figure 82. Combined Thermochemical and Particle Erosion Depths for the MNASA 

Nozzle 
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Figure 83. Impacting Particle Mass Fluxes for the MNASA Nozzle 






Table 10. Erosion and TAD Summary for the MNASA Nozzle 
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4.3 ASRM NOZZLE EROSION 


Erosion depths calculated for the ASRM nozzle are presented in Figure 85. This 
figure presents erosion depth as a function of station. The figure also includes the material 
thickness at each analysis location. Results are presented for both thermochemical 
erosion and combined thermochemical and particle erosion. Station 6 is the nozzle throat. 
Following are comments relative to these results. 

1. The magnitude and characteristics of the thermochemical erosion are consistent with 
those predicted and measured in the MNASA nozzle. 

2. The augmented heating at the exit cone material interface does not appear to cause 
excessive erosion. 

3. The particle erosion in the throat entrance regions is extremely severe and probably 
excessive. 

4. The particle erosion in the nozzle exit cone is not calculated to be excessive. 

This erosion was calculated based on the 'backed out" mass fluxes from the MNASA 
nozzle. Figure 86 compares the shape for the two nozzles. This figure presents surface 
angle as a function of normalized axial position. The surface angle is presented beginning 
at the nozzle throat and continuing to the exit plane of the nozzle. The axial position is 
normalized based on the axial length from the throat to the exit plane. Since the nozzle 
contours are not identical, the mass flux results were applied to the ASRM nozzle based 
on the local surface angle, i.e., an MNASA exit cone mass flux level was applied at the 
location in the ASRM nozzle which had the same surface angle. This figure shows that 
the ASRM would be expected to experience less severe erosion than the MNASA and it 
would occur further aft. 

Figure 87 presents the TAD's calculated for the ASRM nozzle. Again results are 
presented at the end of the motor firing and after two minutes of cool down. This plot 
also includes the local material thickness for comparison. A summary of the calculated 
erosion and TAD results is presented in Table 11. This table also includes the safety 
margin calculated for each of the analysis positions. For these calculations the safety 
margin is defined as follows: 


safety margin = 


n }-i 

2xe + 1.25xc + di. 


where, 

t; = liner thickness 

e = erosion depth 

c = char depth at the end of bum 

dj = depth to the 100° F isotherm for analysis stations 1 through 13 and 

depth to the 350° F isotherm for analysis stations 14 through 17. 

Based on this definition of safety margin, it appears the design is inadequate throughout 
the nozzle. This is especially true in the entrance regions of the nozzle where particle 
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impact effects were calculated to be the most severe. For the majority of the remaining 
analysis positions the safety margin is only a small negative number (or positive in some 
cases) and is therefore not as great of a concern. 

4.4 ASRM FLEX SEAL RESPONSE 

The response of the flex seal for the ASRM nozzle was analyzed using the 
polyisoprene model described in Appendix B and the convective heating conditions 
described in section 3. As described in Appendix B, the two dimensional conduction 
effects of the polyisoprene/carbon phenolic configuration do not appear significant. 
Therefore, these material response calculations were performed with the CMA computer 
code rather than with ASTHMA. Results from this calculation are presented in Figure 88. 
This figure presents the recession response for both the polysioprene and the carbon 
phenolic. The polyisoprene is calculated to recede at very rapid rates. In fact, the 
calculation shows the polyisoprene should ablate completely through in about 80 seconds. 
This calculation is probably conservative in several respects: 

1. The polyisoprene model was developed from the MNASA motor data for which 
there was probably particle erosion present. Particle impact conditions for the 
ASRM flex seal location should be significantly less. Because of this, the model may 
be somewhat conservative and tend to over predict recession in a pure 
convective/radiative heating environment. 

2. The modeling of the polyisoprene does not account for the substantial differential 
recession which will occur between the polyisoprene and the carbon phenolic. This 
difference in recession can allow the carbon phenolic to shield the polyisoprene from 
both the convective and radiative heating environments. This effect was not taken 
into account in these analyses. 

In spite of these conservative aspects of the calculation, it is obvious that the flex seal is a 
potential problem area and needs additional attention. Section 6 provides some 
recommendations for addressing the flex seal. 
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Figure 85. Calculated Erosion Depths for the ASRM Nozzle 
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Figure 87. TAD for ASRM Nozzle 
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Figure 88. ASRM Nozzle Flex Seal Recession 
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Table 11. Erosion, TAD and Margin of Safety Summary for the ASRM Nozzle 



5.0 CONCLUSIONS 


In conclusion, this study has been extremely beneficial in increasing the 
understanding of the MNASA motor and applying its results to the development of the 
ASRM motor. Several key aspects of the MNASA calculations include: 

1. The analytical procedures used in performing these analysis do an excellent job of 
predicting the heating and material response for the MNASA nozzle. 

2. The method developed to account for the augmented heating at the exit cone 
material interface does an adequate job of predicting the measured material response. 

3. The particle flow calculations and mechanical erosion predictions used in this study 
represent an improvement over standard procedures used in the solid rocket motor 
industry. Only the larger particle sizes (120 pm) are calculated to impact the nozzle 
in the throat entrance region. The calculated results were found to be consistent with 
the measured data. 

4. No particle impacts were calculated to occur in the nozzle exit cone. However, the 
measured data suggests impacts do occur. To account for this modeling deficiency 
an empirical approach was used to "backout" impacting particle mass fluxes. 

5. The TAD’s calculated for the MNASA motor compare very well to the measured 
values for the MX4996 material. The TAD's are significantly underpredicted for the 
FM5939 material. 

6. A material response model was developed for the polyisoprene flex seal material 
from the MNASA motor firing data. This model treats the polyisoprene as a surface 
sublimer. However, because of uncertainties in the MNASA environment, this model 
may be somewhat conservative when applied to the ASRM motor. 

Using the models developed and the lessons learned from the MNASA nozzle 
analysis, similar analyses were conducted for the ASRM nozzle. Based on these 
calculations the following conclusions can be drawn: 

1. Particle erosion will cause excessive material removal in the throat entrance region 
and therefore, large negative margins of safety. 

2. Augmented heating at the exit cone material interface will cause only slightly 
excessive erosion (based on small negative safety margins). 

3. Particle impacts in the exit cone will not cause excessive erosion (safety margins of 
about zero). 

4. The polyisoprene in the flex seal is calculated to ablate completely through in only 
80 seconds and appears to be a significant problem area. 
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6.0 RECOMMENDATIONS 


There are several aspects of the results which suggest that further studies are 
needed to better quantify margins of safety for the ASRM nozzle. The following 
recommendations, in order of importance, describe areas which require additional study. 

1 . Additional effort needs to be directed toward understanding the polyisoprene 
material response. This should include testing of the polyisoprene in controlled 
environments so that the material recession mechanisms can be better understood. 

At this point it is not clear if the material behaves as a surface sublimer or if the shear 
environment plays an important role in the material removal. Additionally, the role 
of thermochemical erosion is not known. The best place to quantify these effects is 
in arc plasma generator tests where the convective heating and thermochemical 
environments can be controlled and made to duplicate the environments of the 
ASRM nozzle. 

2. The particle erosion of the throat entrance region needs to be better understood. In 
order to further validate the procedures used in this study, the same procedures 
should be evaluated against the RSRM motor firing data. This motor in many 
respects does a much better job of duplicating the erosion environment of the 
ASRM than does the MNASA motor. Comparing to data from this motor should 
add a great deal to understanding the erosion in the ASRM nozzle as well as 
understanding the accuracy of the analysis procedures. 

3. Material response models for the MX4996 and FM5939 need to be developed so that 
margins of safety, especially with regard to char depth can be more accurately 
predicted. These models should be developed with the use of analog and arc heater 
testing so that material response data from representative and well quantified 
environments are used to understand the material response. This model development 
is especially needed for the FM5939 since there was such a large difference in 
calculated and measured TAD's for this material. 
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APPENDIX A 


INTERIM THERMAL RESPONSE MODELS FOR MX4996 AND FM5939 

CARBON PHENOLICS 

This Appendix presents the interim thermal response models which have been 
developed to model the standard density MX4996 carbon phenolic and the low density 
FM5939 carbon phenolic used in the MNASA ASRM 48-5 nozzle and the full scale 
ASRM nozzle. These models used thermal conductivity and TGA data generated at 
Southern Research Institute (SRI), Birmingham, Alabama. Specific heat of the materials 
was estimated based on similarity to other carbon phenolic materials. Additionally, 
chemical composition of the virgin materials was measured under this program. These 
measurements were made at Galbraith Laboratories in Knoxville, Tennessee. The 
compostion of the pyrolysis gases have been inferred from these measurements of virgin 
material composition and the assumption that the char materials are composed entirely of 
carbon. It should be stressed that these are only interim models and that more data needs 
to be gathered before thermal response calculations can be made with a high level of 
confidence. This data should be in the form of analog tests at Southern Research Institute 
and arc plasma tests conducted at Aerotherm. 

The interim thermal response model developed for MX4996 carbon phenolic is 
presented in Tables A-l through A-6 and Figures A-l through A-3. Data is presented in 
both tabular and graphical form. The char thermal conductivities used in this model were 
estimated based on dynamic thermal conductivities determined by Aerotherm for other 
similar carbon phenolics. These estimated values were used in place of the laboratroy 
measured values since the laboratory values traditionally underpredict thermal response 
data by a significant amount. A comparison of the estimated and measured values is 
provided in Figure A-3. 

The interim thermal response model developed for the low density FM5939 carbon 
phenolic is presented in Tables A-7 through A- 12 and Figures A-4 through A-6. The 
conductivities in this case are taken from the SRI measured data. Estimated conductivities 
were not developed since the low density material was believed to be significantly different 
from other materials for which models were available. 
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Table A-l. Decomposition Kinetics for MX4996 Carbon Phenolic 


Reaction;: 

;||l)ensily ' 

Residual 

'' • ' 

Pre-]Expo 

: \ ' 
Xv'.'X'X'i’MvM v 

. Is ■ : ■/ s: ■ 

'"L? “ - 

1 

WZERO(I) 

(Lbm/Ft 3 ) 

RR(I) 

(Lbm/Ft 3 ) 

B(I) 

(1/SEC) 

PSI(I) 

Q 

E(I) 

(DEGR) 

TRES(I) 

(DEGR) 

A 

3.6202 

.00000 

334.19 

1.2900 

9702.2 

540.00 

B 

2.5668 

.00000 

.51330E+12 

1.5200 

39034 

1008.0 

C 

169.45 

146.55 

15.590 

1.2200 

12297 

729.00 


Initial Composite Density WINTL 

GAMA 


= 87.818 (Lbm/Ft 3 ) 
= .500 


TGA Data for MX 4996A 

20 Deg C/Min Heating Rate 


Density (Lbm/Ft A 3) 



Figure A-l. Comparison of Measured and Calculated TGA Data 
for MX4996 Carbon Phenolic 
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Table A-2. Elemental Composition for MX4996 Carbon Phenolic 


■■nhhh 




% Carbon 

88.52 

58.94 

100 

% Nitrogen 

0.50 

3.20 


% Oxygen 

3.65 

23.19 


% Hydrogen 

2.31 

14.67 


% Undetermined 

5.02 




Table A-3. Estimated Specific Heat for MX4996 Carbon Phenolic 


M : - : Temperature W$MM; 

Virgin (Btn/Lbm-R) 

Char (Bt»r/Lbm-R) 

530 

0.26 

0.26 

824 

0.35 

0.31 


0.37 

0.35 

1500 

0.43 

0.41 


0.45 

0.44 

2500 

0.46 

0.47 

3000 

0.48 

0.49 

5460 

0.51 

0.51 



Figure A-2. Estimated Specific Heat for MX4996 
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Table A-4. RESIN DATA 


Type: SC1008 

Mass Fraction: 35% (Assumed) 

Resin Heat of Formation: -1080 Btu/Lbm (Assumed) 
Material Heat of Formation: -378 Btu/Lbm 


Table A-5. Virgin 1 

rhermal Conductivity for MX4 

•996 Carbon Phenolic 

111 Temperature (RF #1 

Across Ply (Btu/Ft-Sec-R) 

Ifill (Btu/Ft^S^c-R) - 1 - 

460 

1.505E-04 

2.373E-04 

610 

1.505E-04 

2.373E-04 

810 

1.505E-04 

2.732E-04 

935 

1.505E-04 

2.917E-04 

1410 

1.505E-04 

2.917E-04 

5460 

1.505E-04 

2.917E-04 


NOTE 1: The values listed above are measured values up to a temperature of 1000 F 


Table A-6. Estimated Char Thermal Conductivity 


Temperature (R) i 

Across Ply Conductivity 

III Fill Conductivity H 
'?• (BtWFt-Sec-R^ ; ^l 

460 

9.64E-04 

2.33E-03 

685 

9.64E-04 

E-2.3303 

985 

1.22E-03 

2.74E-03 

1935 

1.27E-03 

2.63E-03 

2960 

1.27E-03 

2.43E-03 

3410 

1.36E-03 

2.72E-03 

4260 

2.09E-03 

4.67E-03 

5160 

2.77E-03 

6.90E-03 

5460 

3.00E-03 

7.73E-03 
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Estimated Dynamic and Standard Measured Values 



Table A-7. Decomposition Kinetics for FM5939 Low Density Carbon Phenolic 


Reaction 

Initial 

1 ::j Den$ityJ|| 

Density 

Factor 

: I)^%act9r; j | 

Act Energy ; 
' Factor : f| 

||Min;;React:| 

I 

WZERO(I) 

(Lbm/Ft^) 

RR(I) 

(Lbm/Ft^) 

BCD 

(1/SEC) 

PSI(I) 

(-) 

E(D 

(DEGR) 

TRES(I) 

(DEGR) 

A 

2.1890 

.00000 

917.51 

1.3100 

10646. 

540.00 

B 

8.9951 

.00000 

.20191E+10 

2.2700 

32106. 

927.00 

C 

114.74 

100.09 

5026.1 

1.6600 

19562. 

927.00 


Initial Composite Density WINTL= 62.962 (Lbm/Ft^) 

GAMA = .500 


A5 



TGA Data for FM5939 LDC 
20 Deg CfMin Heating Rate 



Temperature (K) 

Figure A-4. Comparison of Measured and Calculated TGA Data for 
FM5939 Low Density Carbon Phenolic 


Table A-8. Elemental Composition for FM5939 Low Density Carbon Phenolic 


Element Ifll 

: ^ Virgin 

l ; n :p$ Siiiiiiiiiiii 

Char 

% Carbon 

86.39 

55.55 

100 


0.83 

4.23 


% Oxygen 

5.01 

25.69 


% Hydrogen 

2.83 

14.53 


% Undetermined 

4.94 




Table A-9. Estimated Specific Heat for FM5939 Low Density Carbon Phenolic 


Temperature (R) ■ mm 

Virgin 

(Btfo&bm-R) 

lllilCfia 

ii(Btu/Lbm-R> 

530 

0.26 

0.26 

824 

0.35 

0.31 

1000 

0.37 

0.35 


0.43 

0.41 


0.45 

0.44 

2500 

0.46 

0.47 


0.48 

0.49 

5460 

0.51 

0.51 
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Figure A-5. Estimated Specific Heat of FM5939 Low Density Carbon Phenolic 


Table A- 10 Resin Data 
Type: 91 LD 

Mass Fraction: 35% (Assumed) 

Resin Heat of Formation: -1080 Btu/Lbm (Assumed) 
Material Heat of Formation: -378 Btu/Lbm 


Table A-l 1. Virgin Thermal Conductivity for FM5939 Low Density Carbon Phenolic 


M 

• Temperature 1 (R) . ||f§| 


Fill Coudtodlvity 

■ illiaiiliiiil 

460 

8.102E-05 

1.331E-04 

590 

8.102E-05 

1.331E-04 

760 

9.722E-05 

1.528E-04 

960 

9.954E-05 

1.713E-04 

1260 

1.065E-04 

1.852E-04 

1460 

1.100E-04 

1.910E-04 

5460 

1.100E-04 

1.910E-04 


NOTE 1: The values listed above are measured values up to a temperature of 1000 F 
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Table A- 12. Char Thermal Conductivity for FM5939 Low Density Carbon Phenolic 


Temperature (R) 


v.Acros; 

mi# 


Missal 

460 

♦5.505E-05 

1.331E-04 

630 

♦21.063E-05 

5.093E-04 

1200 

♦30.862E-05 

6.944E-04 

1685 

♦34.103E-05 

7.639E-04 

2610 

♦35.610E-05 

6.481E-04 

3460 

♦34.376E-05 

6.944E-04 


♦38.732E-05 

8.333E-04 

5260 

♦57.422E-05 

14.815E-04 


♦-Estimated value 


NOTE 2: Fill conductivity values were taken from measured data pre-charred at 3500 F 



Temperature (R) 

Figure A-6. Thermal Conductivity for FM5939 Low Density Carbon Phenolic Char 
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APPENDIX B 


POLYISOPRENE ABLATION MODEL 

A material response model has been developed to predict the behavior of the 
polyisoprene flex seal material. An initial model was developed which utilized an 
Arrhenius fit of polyisoprene TGA data. This model treated the material as a 
decomposing ablator. Thermochemistry tables were generated with the ACE module of 
the NAT code. These tables treated the char material as pure carbon and the pyrolysis gas 
as an equilibrium mixture with five moles of carbon to eight moles of hydrogen. To obtain 
a reasonable match of the measured erosion data, a surface fail temperature of 4500° R 
was used. Because of the high fail temperature, this was not believed to be a good 
material model in that excessively high surface temperatures would be calculated. 

Because of the lack of realism of the fail temperature model, this approach was 
abandoned. A new model was developed which treated the material as a surface sublimer. 
With this approach the surface recession rate of the material is assumed to be a function of 
the surface temperature. This behavior was model with a zeroeth order Arrhenius rate 
expression of the form: 


s = Aexp(-B/T w ) 


where, 

s is the surface recession rate 

A is the pre-exponential factor 

B is the activation energy 

T w is the wall temperature 

For this model ,the material was assumed to ablate between temperatures of 1 500° R and 
1600° R. This temperature range was selected based on TGA results which showed the 
material to be fully decomposed at approximately these temperatures. Using this approach 
the following Arrhenius expression was developed: 

s = 1.0037E+29exp(-110,530/T w ) 

With this expression, pseudo thermochemistry tables were developed which are 
compatible with the CMA computer code. Additional details on this procedure are 
contained in Reference B-l. CMA still uses a surface energy balance technique to 
calculate recession rate; however, it is significantly simplified compared to the most 
general case. The energy balance which is performed is illustrated in Figure B-l. This can 
be written in equation form as : 

Qconv + Qrad - (pS) w Hw - 06T W 4 + m c H c = qoond 
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All of these terms are defined as discussed in the CMA User’s Manual with the exception 
of H w which is given by: 


Hw = 


Heas + B'Hi. 


__ J-igas 


l + B' 


where, 

Hgas = enthalpy of propellant edge gases 
Hj so = enthalpy of polyisoprene gases 

A separate routine has been written which generates the necessary B’ tables. Instead of 
non dimensional pyrolysis gas rate B g ’ being an independent variable, the transfer rate 
h/Cp is the independent variable. The last important aspect in the model development is 

the determination/selection of a heat of formation, AH^ . The heat of formation is used in 
determining the enthalpy of the char material, H c is given by the following 
expression: 


H c = AH° f +J Tw C I 4T 

With this approach, a heat of formation was determined in an iterative fashion. A value of 
-4500 BTU/lbm was found to give reasonably good agreement with the polyisoprene 
erosion data taken from the MNASA ASRM 48-5 motor firing. Results using this model 
are presented in Figures B-2 and B-3. Figure B-2 presents the calculated erosion as a 
function of time. Total erosion is calculated to be slightly greater than 0.95 inches. This 
is approximately the amount measured in the ASRM 48-5 firing. Figure B-3 shows 
calculated temperature profiles within the material at various points in time. Notice the 
extremely thin heat affected depth which is calculated. This is because of the conductivity 
levels selected to model this material. Figure B-4 shows conductivity as a function of 
temperature. These values are for an ASRM internal insulator (RP7136) and were taken 
from Reference B-2. Of course, conductivity values above about 1600° R are meaning 
less in that these material temperatures are never attained. 

Two dimensional conduction effects were considered a possible contributor to the 
polyisoprene material response. To investigate their effects, a two dimensional ASTHMA 
model of the flex seal specimen for the MNASA motor was constructed. This model is 
presented in Figure B-5. This figure shows the grid used to perform the calculation. The 
model includes both polyisoprene and carbon phenolic materials. These materials are 
contained in alternating columns of the grid. Thermal and recession calculations were 
performed with this grid to calculate the polyisoprene results. This calculation yielded 
0.969 inch of recession for the polyisoprene. This amount is only 1.4% greater than the 
0.956 inch calculated with the one dimensional approach. Based on this result it was 





decided that two dimensional effects were insignificant and therefore all remaining 
calculations for the polyisoprene were performed with the one dimensional CMA code. 

It should be stressed that this model is a "pure fabrication" in that it is only a best 
guess at this point. Hopefully, this best guess model will have proper sensitivity to the 
important parameters such that it can be used with reasonable accuracy for the full scale 
motor. Again, it should be stressed that this model is merely an interim and in no way 
should it replace the material characterization effort which is needed to better understand 
and model this material's behavior. 
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Figure B-l. Polyisoprene Surface Energy Balance 
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